Supercurrent-induced Peltier-like effect in superconductor/normal- metal weak links 
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The local nonequilibrium quasiparticle distribution function in a normal-metal wire depends on the 
applied voltage over the wire and the type and strength of different scattering mechanisms. We show 
that in a setup with superconducting reservoirs, in which the supercurrent and the dissipative current 
flow (anti)parallel, the distribution function can also be tuned by applying a supercurrent between 
the contacts. Unlike the usual control by voltage or temperature, this leads to a Peltier-like effect: 
the supercurrent converts an externally applied voltage into a difference in the effective temperature 
between two parts of the system maintained at the same potential. We suggest an experimental 
setup for probing this phenomenon and mapping out the controlled distribution function. 
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Many of the well- understood phenomena in mesoscopic 
physics can be probed within the linear response of 
a physical system to an applied external perturbation, 
i.e., they are governed by equilibrium physics. Recently 
the attention has turned more towards the study of ef- 
fects far from equilibrium. The quasiparticle distribu- 
tion function f(x;E) characterizing the nonequilibrium 
was measured in a normal-metal (N) wire between two 
large reservoirs^ through a superconducting (S) tunnel 
probe. This yielded useful information on the residual in- 
teractions between the Fermi-liquid quasiparticles. This 
nonequilibrium distribution was used to control the su- 
percurrent in a normal-metal weak link^*^ Both of these 
setups serve as different types of local probes for f(x; E). 

As a further step, we describe the control of f(x;E) 
via the supercurrent. We show that, unlike other control 
parameters, it changes the profile of the effective temper- 
ature through the sample in the form of a large Peltier 
effect, i.e., heating the electrons in one part of the struc- 
ture, and cooling them in another — even in the case 
of complete electron-hole symmetry. Moreover, we show 
how the two types of measurements for f(x; E) can be 
combined within the same sample. 

We concentrate on studying a diffusive normal-metal 
wire, where elastic scattering is the dominant scattering 
mechanism. In the absence of superconductivity and for 
wires much shorter than the inelastic scattering length, 
the steady-state distribution function between two reser- 
voirs with chemical potentials \i\ and p,i has a double- 
step form, interpolating between the two Fermi functions 
in the reservoirsi 

When the N reservoirs are replaced by superconduct- 
ing ones, the leading transport mechanism at energies 
below the superconducting gap A is Andreev reflection. 6 
This leads to a penetration of superconducting corre- 
lations into the N wire (superconducting proximity ef- 
fect). It modifies the charge and energy conductivities 
and we may introduce the corresponding diffusion coef- 
ficients Vt(x;E) and T>l(x;E) depending on space and 
energy^ More importantly, the proximity effect allows su- 
percurrents to flow through the N wire. To describe these 
effects, it is convenient to separate fix; E) into symmet- 



ric and antisymmetric parts relative to the chemical po- 
tential us of the superconductor, 

f T (E) = l-f(ps-E)-f(n s + E), (1) 
f L (E) = f( f x s -E)-f(ns + E). (2) 

Below, we choose us = 0. These functions describe 
charge and energy distributions, respectively. They sat- 
isfy the kinetic equations^ 

^7 = 0, j T = -D T (x)^+j E f L + T(x)d x f L] (3) 



T(x)d x f T . (4) 



Here we assume no energy relaxation, so the kinetic equa- 
tions describe the conservation of jt{E) and Ejl, the 
spectral charge and energy currents, respectively. Terms 
T>t, "Dl, je, and T can be found from quasiclassical equa- 
tions for the retarded Green's function in the diffusive 
limitiS All of them depend on the phase difference <j> 
between the superconductors such that for = 0, je 
and T vanish. In our case, the charge diffusion coeffi- 
cient T>t is increased at most up to 20% from its normal- 
state value T>t = lr^ whereas for energies below A, Dj, 
tends towards zero near the S interface, effectively pro- 
hibiting energy transport. The term T{x; E, (f) (Ref. @) 
is obtained as a cross term from the retarded and ad- 
vanced Green's functions. In general, it is much smaller 
than the other coefficients. The supercurrent is de- 
scribed by a spectrum je (E; <fr) of supercurrent-carrying 
states ; 11 ' 12 ' 13 which yields a contribution j_e/l(cc) to the 
spectral charge current and, under nonequilibrium con- 
ditions involving Jt{x) ^ 0, a contribution to the energy 
current Ej E fT(x). 

These kinetic equations have to be supplied with 
boundary conditions. At N reservoirs, electrons are sim- 
ply transmitted and the distributions have to match 
Fermi functions with shifted chemical potentials. At the 
NS interface for \E\ < A, Andreev reflection prohibits 
the transfer of energy into S yielding jr, = 0. The charge 
distribution is continuous, which leads to friE) = at 
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the NS interface assuming that there is no charge imbal- 
ance in the leads. 

The nonequilibrium distribution function may be char- 
acterized through its moments, the local chemical poten- 
tial fJ>(x) and the local effective temperature T e ff(x). The 
previous characterizes the charge distribution function as 
fx(x) = J °° dE fx(x;E). The effective temperature de- 
scribes the amount of heat in the electron system and is 
related to the energy distribution function via 

2 r poo 

-^T 2 s (x) = J dEE[f L , (x; E) - f L (x; E)}, (5) 

where £ — (ir 2 kg/3e 2 ) is the Lorenz number and the 
corresponding zero-temperature distribution has a step- 
function form fL,o(x; E > 0) = 1 — 6[E — fi(x)]. 

In the absence of the supercurrent, the kinetic Eqs. 
and (@J are not coupled and, consequently there is 
no thermoelectric coupling between the applied voltage 
and the energy currents. This results from the assump- 
tion of bands with complete electron-hole symmetry in 
the derivation of the formalism. Beyond the limits of 
the formalism, it is known that electron-hole symmetry 
breaking leads to small thermoelectric effects in normal 
metals, limited by the tiny factor ksT/ep^ 

Below, we study a multiterminal setup depicted in 
the inset of Fig. varying the voltage between the 
N and S reservoirs while maintaining the superconduc- 
tors at equal potentials allows one to vary the distribu- 
tion function in the phase-coherent wire. Such a device 
has already been implemented for controlling the criti- 
cal current for the dc Josephson effect. 5 It permits to 
study the supercurrent under nonequilibrium conditions 
without the complications caused by the ac Josephson 
effect and is, hence, an appropriate system for demon- 
strating the physics outlined above: As the energy flow 
EjEfrix) carried by the extra quasiparticles injected 
into the supercurrent-carrying states cannot pass into the 
superconductors, it has to be counterbalanced by another 
energy flow. This flow is driven by the gradient of the 
energy distribution function EV^OxfL and hence, the ap- 
plied control voltage is converted into a gradient of the 
effective temperature through the supercurrent. 

Solving Eqs. © and @ for = and E < A is 
similar to a two-probe N-S case£: Jl stays constant 
throughout the phase-coherent wire at its value in the 
N reservoir, f° L {V) = (tanh[(£ + eV)/2k B T] + tanh[(E - 
eF)/2fceT])/2 and fx is slightly modified from the linear 
space dependence due to the proximity effect on T>t 10 
Increasing the phase <fi induces a finite supercurrent into 
the weak link, thereby coupling and /t- First neglect- 
ing the small coefficient T, we get 

<% h d / df T \ 2 h 

dx-{ VT -dx-) =jE v2- (6) 

Assuming that je is small, we observe that the major 
change due to the supercurrent is expected for x); 
particularly, it will depend on space. 




FIG. 1: Quasiparticle distribution function f(x;E) in the 
right-horizontal arm for voltage V = 20-Bt/c, temperature 
T = AEx/kB, and phase difference <f> = n/2 between the su- 
perconductors. The large deviations from the rounded stair- 
case form are created by the supercurrent flowing in the struc- 
ture. 

In general, a closed-form solution for Jl(x;E), 
Jt{x]E) cannot be found. Therefore, we focus on a nu- 
merical solution of both the spectral and kinetic equa- 
tions. Here and below, we assume that all the ener- 
gies are below A. The effect of the supercurrent on the 
distribution functions is clearest at a low temperature 
ksT < Et- The resulting distribution function f(x,E) 
for the right-hand horizontal arm is plotted in Fig. ^ f° r 
(j) = 7r/2, yielding a supercurrent close to its maximum. 
As expected, the antisymmetric part of f(x] E) has be- 
come space dependent, its energy dependence following 
that of je- Fixing a position in space, chosen, for exam- 
ple, near the NS interface in the left-hand side horizontal 
arm, allows us to observe how the distribution function 
changes as a function of phase 0, i.e., as it is driven by 
the supercurrent. This is illustrated in Fig. [21 where 
Jl(E; 4>) is plotted for a few values of <\>. 

In the three-probe case, the chemical potential n{x) 
interpolates nearly linearily between the chemical poten- 
tials of the superconductor and the normal reservoir and 
varies only little with the supercurrent. The changes in 
the effective temperature T e ff are much more pronounced. 
In the absence of the supercurrent, T c q is 

T c ° ff ee T cff (0 = 0) = + {V 2 (M*)/e) 2 }/£o. (7) 

Both T c g(x; 4> = 0) and [i(x; <fi = 0) are symmetric in the 
two horizontal arms. The supercurrent-induced change 
in Jl{x;E) can be described through the change of the 
effective temperature compared to Eq. {7J), such that 
T oS (x) = [T° s (x) 2 + S(x; V) + 5fi(x; cj))] 1 ' 2 , where 

(\ f'°° 

S(x;V) = — T dEE(f L (E;V)-f L (x;E)) (8) 

7T K B J Q 

and Sfi(x; (j>) = [n(x; <p) 2 — fi(x; (f> — 0) 2 ]/2 describes the 
change in the local chemical potential due to the super- 
current, a much smaller term than S(x). The kinetic 
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FIG. 2: Supercurrent-driven distribution function /l(-E) at 
the left NS interface as a function of energy for cf> — (solid), 
4> = 0.12tt (dotted), (f> = 0.24vr (dashed), and <j> = tt/2 
(dash dotted). The result is obtained with T = AEr/ks 
and V = 20Er/e. The corresponding changes of /t by the 
supercurrent in the right arm have the opposite sign. Inset: 
the system under consideration. We assume symmetric hori- 
zontal wires of length L/2. This length defines the Thouless 
energy of the weak link, Et = TlD/L 2 . The resistance of the 
weak link is i? w and of the vertical wire R c . Measurement 
of the predicted effects can be performed by placing a su- 
perconducting tunnel probe at position x — x c , near the NS 
interface. 




FIG. 3: Supercurrent-induced change S(x; V) in the effective 
temperature as a function of voltage eV/Ex at different po- 
sitions in the weak link with <j> — tt/2. From top to bottom: 
x — 0, x = L/4, x = L/2, x = 3I//4, and x = L. Here x — 
and x = L correspond to the left and right S interfaces and 
x — L/2 to the crossing point. Inset shows the phase de- 
pendence of S(x = 0) for eV = 12E T (solid) and eV = 8E T 
(dashed). In both curves, the bath temperature T = 0. 



equations imply that the supercurrent-induced change of 
the distribution function //, is antisymmetric between 
the two arms, hence so is S(x), i.e., T e ff increases in one 
arm and decreases in the other one. Hence, the system 
works analogously to a Peltier device, where the control 
current is replaced by the supercurrent: the supercur- 
rent "cools" one part of the system, transferring the heat 
to another part. The supercurrent-induced temperature 



change S(x; V) is illustrated in Fig. |3| 
To obtain an estimate for S(x; V, <p) 
V T (E;x) = 1 and find 



we approximate 



2E^R V 



L2(R W + R C ) 



dEEj E (E)f°(E) 



L/2 



'daf 



V L (x';E) 
(9) 

At low temperatures ksT <C eV, fr reduces to a step 
function around the potential eV in the reservoir, cutting 
the integration off at E — eV. Thus, this current-induced 
temperature change, which is similar to the Peltier effect, 
is much larger than in conventional single-metal setups. 

For the measurement of the predicted effects in the 
distribution function, we suggest a setup shown in the 
inset of Fig. |2 - very similar to those used in Refs. Ill'2t 
Such a setup has also been used as a local thermometer 15 
of the electronic temperatures. A superconducting wire 
is connected to the horizontal arm via a highly resistive 
tunneling layer (I) at position, say, x = x c . The dc cur- 
rent is then given by the tunneling quasiparticle current 



7,7 eR 



dEps (E + p)p(E) [fo (E+fx)-f(E)], (10) 



where N s ps(E) is the BCS density of states (DOS) of 
the tunnel probe, N^p(E) is the local DOS in the meso- 
scopic wire at x — x c , N$ and iV/v are the normal- state 
DOS's for the two materials at E = Ep, fo(E) is the 
Fermi function, and f(E) is the distribution function to 
be measured. When all the wires are in the normal state, 
the resistance through the tunnel junction is R. We can 
separate this expression as Ij = I\ + I2, where I\ is 
the tunneling current for the equilibrium system V = 0, 
probing p(E) and 

1 r°° 

h = — / dEp(E){f T (E)[p s (p + E)+ Ps (p - E)] 
ere j Q 

+ {f L (E) - tanh(£/2fc B T)] \p s {fM + E)-p s (p- E)}} 

(11) 

depends on the state of the wire, and for an equilibrium 
state, V = 0, vanishes. In order to isolate I2, one can first 
determine I\ as a function of the supercurrent by investi- 
gating the equilibrium case. Then, I\ may be substracted 
from the nonequilibrium results, leaving only currents I 2 . 
Moreover, 7a (m) + h{— m) is proportional to the first part 
of Eq. lfTT|) . dependent on fr(E), and hifj) — p) to 
the second, dependent on Jl{E). 

With the above setup, the distribution functions may 
be characterized as a function of both the voltage V 
and the supercurrent driven through the weak link. The 
setup also makes it possible to measure the local distribu- 
tion function both through the NIS contact and through 
the SNS critical current. These two independent probes 
should permit to distinguish the contributions from dif- 
ferent inelastic scattering effects along the lines of Ref. Q. 

So far we have completely neglected inelastic scatter- 
ing in the wires. We can include energy relaxation due 
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to electron-electron scatterin g p henomenologically, gen- 
eralizing the method of Ref. 16 to include the effect of 
supercurrent. In the limit L > ( £ , we may describe the 
nonequilibrium distribution functions by Fermi functions 
with local chemical potential and temperature. In this 
case, assuming for simplicity T>t = 1 and T = 0, we 
can integrate the two kinetic equations over energy, ob- 
tain kinetic equations for fx(x) and T e g(x) and find in the 
limit of high A 

dln{x) = -dj s {x), (12) 
e 2 X f(x)d x T eS (x) = -p,(x)d x n(x) + Q S - (13) 

Here I s (x) = [J dEj E (E)f L (E,x)}/2 is the local su- 
percurrent, 2C e 2 f = - J dEET> L dTfL and 2fl(x) — 
— J dEEULd^fL describe the local temperature and 
chemical potential modified by T>l, respectively, and 
Qs = [J dEEj E (E)f T (E, x)]/2 is the energy current car- 
ried by the nonequilibrium supercurrent. The first equa- 
tion states the conservation of the total current whereas 
the latter describes the temperature profile. In the ab- 
sence of the proximity effect, these yield the effective 
temperature given in Eq. JJJ). Similarly as above, the 
effective temperature can also in this case be tuned via 
the supercurrent, through the control of Qs- 

The predicted effect resembles a previously studied 
phenomenon in bulk superconductors^ where a tem- 
perature gradient along with a supercurrent generates a 
charge imbalance in S. Here, the finite voltage (described 



through fx) along with the supercurrent produces a tem- 
perature gradient (spatial variation of /l)^ 

In Ref. flfl another thermoelectric effect, the ther- 
mopower, has been measured experimentally in a similar 
type of a system. The coupling of the distribution func- 
tions through the supercurrent may explain part of the 
observed effects. In Ref.[2£j, thermopower has been stud- 
ied in the regime of high tunnel barriers and within lin- 
ear response, leading also to an unexpectedly large effect. 
In that paper, all the distribution functions are, besides 
minor corrections, in quasiequilibrium: the transport is 
essentially driven by the discontinuities at the tunneling 
barriers. Moreover, Ref. |^] studies the Andreev inter- 
ferometers through a numerical scattering approach, and 
predicts an oscillating thermopower as a function of the 
phase. However, there the quasiparticle current and su- 
percurrent do not flow in parallel and the magnitude of 
the effect may be strongly affected by the very small size 
of the studied structure. 

Summarizing, we predict that in a nonequilibrium sit- 
uation created by applying a voltage between a normal 
metal and two superconductors, the nonequilibrium dis- 
tribution functions in the normal-metal wire can be tuned 
by the supercurrent flowing between the superconduc- 
tors. This results in a supercurrent-controlled Peltier ef- 
fect. The predicted effect can be observed by the mea- 
surement of the tunneling current from an additional su- 
perconductor. 
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